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Abstract. In this paper, we study various dissipative mechanics asso- 
ciated with the Boussinesq systems which model two-dimensional small 
amplitude long wavelength water waves. We will show that the decay 
rate for the damped one-directional model equations, such as the KdV 
and BBM equations, holds for some of the damped Boussinesq systems 
which model two-directional waves. 



1. INTRODUCTION 

Considered here are waves on the surface of an inviscid fluid in a flat 
channel. When one is interested in the propagation of one-directional irro- 
tational small amplitude long waves, it is classical to model the waves by 
the well-known KdV (Korteweg-de Vries) equation (see 23 ) 

u t + u x + u xxx + uu x = 0, 

or its regularized version, the so-called regularized long wave equation or 
BBM (Benjamin-Bona-Mahony) equation, 

u t + u x - u txx + uu x = 0. 

When one is dealing with two-directional waves, and the effects of wave inter- 
actions and/or wave reflections are not excluded from the study, a restricted 
four-parameter family of systems (see jH]), 

(1.1) * + «, + M. + «ti«-^ = 0, 

Ut + Vx + uu x + cq xxx - du xxt = 0, 

may be used. The dimensionless variables rj(x, t), u(x, t), x, and t are scaled 
by the length scale ho and time scale (ho/g)? where ho denotes the still 
water depth and g denotes the acceleration of gravity. The variable rj(x,t) 
is the non-dimensional deviation of the water surface from its undisturbed 
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position and u(x, t) is the non-dimensional horizontal velocity at a height 
above the bottom of the channel corresponding to Oho with < < 1. The 
constants a, b, c, d are called dispersive constants which satisfy the physical 
relevant constraints 

(CO) a + b + c + d = \ and c + d = ~(1 - 6 2 ) > 0. 

This class of systems contains some of the well-known systems, such as 
the classical Boussinesq system (o = b = c = 0, d = 1/3) (see for example 
HHEIHI ED 113) and the Bona-Smith system (a = 0, b = d > 0, c < 0) 0. 
It is shown in [6: that a physically relevant system in (Jl.lj) is linearly well 
posed in certain natural Sobolev spaces if the constants a, b, c, d satisfy 

(CI) b > 0, d > 0, a < 0, c < 0, 

or 

(C2) 6 > 0, d > 0, a = c> 0. 

It is also shown in |H[7j that above systems have the capacity of capture the 
main characteristics of the flow in an idea fluid. But when the damping effect 
is comparable with the effects of nonlinearity and/or dispersion, as occurs in 
the real laboratory-scale experiments and in the fields (see [71 1161 IT21 117j ). 
it should be considered in order for the model and its numerical results 
to correspond in detail with the experiments. The full system would be 
the Navier-Stokes equations with a free boundary, which is very difficult 
to handle both theoretically and numerically (cf. |21| IH]). Therefore, it is 
useful to construct simpler model systems which are capable of capturing 
the main properties of water waves under various special circumstances. 

For example, many researchers have studied the dissipative one-way prop- 
agation model equations, such as the dissipative KdV and dissipative reg- 
ularized long- wave equations and their generalizations. As a model to our 
study, we recall the results from |2j for the dissipative BBM equation, 

Ut - Uxxt - vu xx + uu x = 0, 
u(x, 0) = uq(x) 

where v is a positive constant. 

Theorem 1.1. Assume uq is in L 1 (1R) flL 2 (K), then there exists a constant 
C such that 

(1.2) |Ki)|| i2 <c{i + t)-w 

Here L X (R) and L 2 (R) are the classical Banach spaces. A similar result 
holds for the corresponding dissipative KdV equation. See also jTJ, ^H] and 
the references therein. 

In this article, we aim to analyze the effect of dissipation on systems (f 1 . If) 
and study the decay rates of solutions (rj, u) toward zero. We will restrict 
our study to the cases where constants a,b,c,d satisfy (C0)-(C1) or (C0)- 
(C2). The goal of this research is to find the appropriate dissipative term 
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(or terms) which will provide the right amount of energy dissipation for all 
wave numbers while keeping the mass conserved. 

In this article, two kinds of dissipations will be considered: 

Complete dissipation: replacing the ^\ in the right-hand side of <|1-1|) by 
the vector ( ) , and 

\^xx J 

Partial dissipation: replace the ( _ ) in the right-hand side of by the 



vector 







We shall first study the decay rates of solutions to the linearized systems 
supplemented with either complete or partial dissipations. These equations 
read 

Vt + u x + au xxx - br/ xxt = ut] x 

(1.3) 



IXXl 

1H ~\~ Tj x -\- CTj xxx du xx t — u XX i 



with v = or v = 1. Systems which satisfy the dichotomy property in the 
Fourier space: 

• decay as t" 1 / 4 for low frequencies (small £); 

• decay as exp(— fit) or exp(— fi£, 2 t) for high frequencies (large £); 

will be identified and studied. It is shown in Section |3] that the dichotomy 
property will lead to the decay rate t" for ||(?7, u)\\ L 2 xH h (h will be specified 
later) . 

For later use, we shall emphasize (and it is easy to check) that the di- 
chotomy property holds true for two fundamentally different equations: the 
linearized BBM-Burgers equation ut — u xx t + u x — u xx = and the linearized 
KdV-Burgers equation ut + u xxx + u x — u xx = 0. If the high frequency part 
of the solution to a system is damped as exp(— (3t), we say that the system 
belongs to the BBM-Burgers class since solutions to linearized BBM-Burgers 
equation feature this property. If the high frequency part of the solution is 
damped as exp(— /?£ 2 t), which is the case for linearized KdV-Burgers equa- 
tion, then we say that the system belongs to the KdV-Burgers class. The 
low frequency parts of the solutions to the linearized BBM-Burgers and 
KdV-Burgers equations behave in a similar fashion. 

The main result in Section is to classify the linearized systems according 
to this property and to prove that for systems which satisfy the dichotomy 
property, a decay rate comparing to ()1.2j) is valid. We shall also present 
some systems where the decay rates can be arbitrarily small, behaving as 
the solution of 

u t -u xxt + u = 0, u(x,0) = u (x), 
where by Fourier transform, 



t 
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and therefore ||n(.,f)||^2 could decay arbitrarily slow. 

In Sectional we extend the linear theory to nonlinear systems and show 
that the decay rate as (|1.2|) is valid for weakly dispersive systems, i.e. sys- 
tems with b > and d > 0, and for some systems in the KdV-Burgers class 
with total dissipation, which include the KdV-KdV system (a = c = \,b = 
d = 0), with small initial data. In Section |S1 the decay rate with respect to 
L°°-norm is presented and in Section 13 spectral method is used on several 
systems to demonstrate that the rates obtained in Section 0] and Section [5] 
are sharp and the constants involved in the bounds are reasonably sized. 

It is worth to note that there are other methods, such as the energy meth- 
ods (like the so-called Schonbek's splitting method applied to the classical 
Boussinesq system |19| in large dimensions), can be used in proving decay 
rate for solutions of these systems. This line of study will be carried else- 
where. We believe that those methods will be helpful especially in the cases 
when b = d so some Hamiltonian is conserved (see [Hj). 

On the other hand, to remove the smallness assumption on the initial 
data, the authors in L 2^ used a kind of Cole-Hopf transformation (that is 
valid for Burgers equation) and were able to control the extra-terms. We 
do not know if this tour de force (feat of skill) is possible for the systems in 
(jl.lj) with dissipation. 

We complete this introduction by introducing some notations. Through- 
out the paper, the standard notation on Sobolev spaces will be used. The 
L P (IR) norm will be denoted as || • \\lp for 1 < p < oo and the H s norm 
will be denoted as || • \\h s - When several variables are involved, we may 
also set L x for L P (R) to specify that we compute the norm with respect to 
the x-variable. The product space X x X will be abbreviated by X and a 
function f = (/i, $2) in X carries the norm 

l|f|U = (||/l||x + ll/2||x) 1 - 

The Euclidean norm of a vector is denoted by 1*1- We will use C and [3 as 
generic positive constants whose values may change with each appearance. 
Fourier transform of a function / is denoted by either / or J-{f). 



2. Notations and Preparations 

2.1. Some notations. Consider u € {0, 1}. As stated before, we plan to 
first estimate the decay rates of solutions to the linear systems 

, s Vt + Ux + au xxx - br\ xxt = uij xx , 

(2.1) 

1H ~r Vx 1 CVxxx du xx t — U XX j 

when t goes to +00. 

Following 0, we introduce the Fourier multipliers 

1 - a e , 1 - c£ 2 
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Since a,b,c,d satisfy (CI) or (C2), uiu? is nonnegative and we denote 

H=(^y 2 and a = {u 1 u 2 ) 1 ' 2 , 
with the conventional notation § = 1- We also denote 

a = —77 and e 



l + b£ 2 ~ l + d£ 2 ' 

Remark 2.1. When a system satisfying (C2) assumption is the subject of 
the study, lu\ and 0J2 do change signs, but u)\u)2 > 0. 

Definition 2.2. Consider a nonnegative function £ — > £(£). The order of 
k (when it exists) is defined as the number m such that 

when |£| — ► +00. The (pseudo-differential) operator k with order m is 
defined by setting 

ku = v iff ku = v. 
Therefore k maps L? x into H x order ^ ( or ffn i n i j{™ order ( K ) j_ 

Since IJ2.1J) is a linear system, it is convenient to use the Fourier transform. 
Let (fj, u) denote the Fourier transform of (n, u) and set Y = (77, w) with 
w = Hu, then (12.11) reads 



(2.2) Y t + AY = 0, 

where 

A ,t\ = ( va isgn(wi)^ 
^' ^sgn(o>2)^o" £ 

is the symbol of the linear (unbounded) operator in 1)2. Since we are 
dealing with a system, .A(£) is a matrix. 

By multiplying Y* on ()2.2|) and taking the real part, 

^4 [(\mo\ 2 +wt,o\ 2 n 



(2.3) 2dtj 

+ u j a(0m,0\ 2 dti + J e(O\w(t,O\ 2 d£ = 0. 
Since a(£) and e(£) are positive, 

(2.4) E(t):= [ \Y(t,0\ 2 dC 

decays towards as t — > 00, where |y| = (|r7| 2 + |t?| 2 ) 1//2 is the Euclidean 
norm on C 2 . 
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2.2. Linear algebra. We recall some facts from linear algebra and then 
apply them to the dissipative systems (j2.1j) . 



Definition 2.3. Let M be a 2 x 2 matrix in the complex space, the norm of 
M is defined by 

insii \ MY \ 
\\M\\ = sup . 

yec 2 \{0} \ Y \ 

Lemma 2.4. Let p(M) denote the spectral radius of a matrix M and tr(M) 
denote the trace of M, then 

||M|| = p{M*M) 1 / 2 < tr(M*M) 1/2 . 

We are now going to bound E(t) (see (12.21) and (J23J) by using the point- 
wise estimate 

(2-5) |F(t,OI<||e-* A || |y (OI- 

Noticing that the matrix A can be written as A = D + U, where D = 

V ^ ^ I represents the dissipation terms and U = ( ■ P % A ^ s S n (^i)^ (7 
ey r ^ \isgn(uj2 )l;o- 

is skew-symmetric. When D and U commute, the behavior of \\e~ \\ with 
respect to £ is characterized by the behaviors of e and a via 

(2.6) l|e -L4 || < e~' min ^"^)' e ^^. 

But when D and U do not commute, more accurate estimate than (|2,6|) can 
be obtained by studying e~ tA in detail. 

We now recall the following lemma (Theorem 9.28 from |13j). 

Lemma 2.5. There exists a unitary matrix Q (i.e. QQ* = Q*Q = L) such 
that 

A = <5*(o x 2 ) Q - 

where Ai and A2 are the eigenvalues of A, ordered by Re(Xi) < Re{\2). 

As a consequence, one can prove (which is given at the end of this sub- 
section) 

Lemma 2.6. There exists C > such that 

(2.7) || exp(-L4)|| < C ( 1 + \z\ min ( t, 1 ) ) exp(-ti?e(Ai)). 



A2 — Ai| 

It is easy to see that Ai and A2 are the roots of the characteristic equation 

(2.8) A 2 - tr(A)A + det(A) = 
where 

(2.9) tr(A) = Ai + A 2 = va + e > 
and 

(2.10) det(yl) = A : A 2 = vae + fa 2 > 0. 
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We now estimate || exp(— tA)\\ by separating the cases A < and A > 
where A is the determinant of 1)2. 8J) . namely 

(2.11) A = tr(A) 2 - 4 det(A) = (e - vol) 2 - 4£ V. 

Lemma 2.7. For any t > and for any ^£t, 

• when A < (perturbation range), 
(2.12) 

||exp(-L4)|| <C(l + tr(A)t)expf-^^ <Cexp^-^^ t 

• when A > (non-perturbation range), 

(2.13) || exp(— <C (l + 2\£\ammU,-j=\}(5tp(-t\i), 
where X\ satisfies 

det(A) , . / . . , 2det(A)\ 

(2.14) — < Ai < min tr(A), ' 



tr(A) - 1 - V t r ( A ) / ' 

Proof. It is worth to note from Lemma 12.51 that 

(2.15) tr(A*A) = |Ai| 2 + |A 2 | 2 + |z| 2 = v 2 a 2 + e 2 + 2fa 2 . 

When A < (perturbation range): matrix A has two conjugate complex 
eigenvalues Ai and A2 with 

|Aa| = |A 2 |, Re(Ai) = Re(A 2 ) = ^^p, |Ai| 2 = |A 2 | 2 = det(A). 

Using (|2T5|l and then lf2~TT)|) - l|2~5jl leads to 
\ z \ 2 = v 2 a 2 + e 2 + 2£V - 2|Ai| 2 

= z^ 2 a 2 + e 2 + 2£V - 2det(4) = (ya - e) 2 < tr(A) 2 . 



Hence (|2.12|) follows from Lemma 127 

When A > (non-perturbation range): tr(A) > and det(A) > imply 
that the matrix A features two real eigenvalues < Ai < A 2 . Then (|2.15|) 
leads to 

| z |2 = v 2 a 2 + £ 2 + 2 ^ a 2 _ tr(A) 2 + 2 = ^2 ff 2 

and H2.13J) is proved by using Lemma l2~o1 Since A 2 < tr(A) < 2A 2 , one sees 
immediately that 

det(A) det(A) 2det(A) 

l^KY ~ l ~^2~~ tr(A) 



and lj2.14j) follows. □ 
It is noted that when A < 0, the dissipation can be considered as a 
perturbation term with respect to the skew symmetric operator. More pre- 
cisely, the decay is the same as pretending U and D commute, up to a linear 
correction. 
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When A > 0, this is no longer valid. In the first case, matrix A has conju- 
gate complex eigenvalues. In the latter case, A has real positive eigenvalues 
and the smallest one monitors the decay estimate. A = is the bifurcation 
point. 

For the sake of completeness, we now give the proof of Lemma 12.61 
Proof of Lemma 12.61 Straightforward computations lead to 

/„-tAi „ e- tA l-e- fA 2 \ 

e~ tA = O* I A i~ A 2 O 

Q { )^ 

where 



Ai — A2 
Lemma 12.41 then yields 

||e- M || 2 <tr(e- tA V tA 



P tXl _ p — t\2 

= - te - tx \ if A 1 = A 2 



(2.16) 

Therefore, if | Ai — A 2 | > 



2tRe(\ 1 ) + e -2tRe(X 2 ) + | z |2 



2 



e — tAi _ g— 1\2 



Ai — A9 



|| e -*A||2 < (2 + 1*1 ) e -2^e(A0 
|Ai — A2P 

which proves part of the lemma. 

Now, for |Ai — A2I > A, where A > will be chosen later, 

\ P — tAl „—t\2 1 9 

(2-17) ^ --^ < ±e- tRe ^; 

|Ai-A 2 | A 

and for |Ai — A2I < A, using |e^ — 1| < \(\ exp(|C|) for any complex number 

c, 

| e — *Ai _ e -U 2 | = e -tRe(Ai)| e -t(A 2 -Ai) _ 1 | 

< e - tRc ( Al )t|A 2 -A 1 |e iA . 



(2.18) 



Therefore, choosing A = ^ in (f2~TT)) and l|2.18|) . 

\„ — t\l _ p — t\2 I 

' I < Cte -tfle(Ax). 



|Ai-A 2 | 

Substituting above into ()2.16j) completes the proof. □ 



3. Decay rate of linear systems 

In subsections l3.1ll3~21 and l3~3l low-frequency (|£| close to 0), high-frequency 
(large |£|) and middle range frequency analysis for the linear systems are per- 
formed respectively. We will identify systems for which there exist positive 
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constants (3 and 5 m such that for any t > and 

• for 1^1 < <5 mj ||exp(-a)|| < C exp(-/^ 2 ), 

' • for |e| > <5 TO , || exp(-U)|| < Cexp(-(3t). 

Here || exp(— tA)\\ is the norm of the linear operator exp(— tA(^)) acting on 
C 2 . The generic constants C and (3 are independent of t and £. If 5 m = +oo 
is feasible in (|3.1|) . the system is in the Kd\ 7 -Burgers class. Otherwise the 
system is in the BBM-Burgers class. A summary of decay rates for the 
linear systems is given in subsection 13.41 

Remark 3.1. It is worth to note that \S. 1]) is sufficient but not necessary 
for proving the desired decay rate. An example is given in Section \S. ,91 where 
\S. 1]) is not valid but the linear system has the desired decay rate. 

3.1. Low frequency analysis. We now prove that for |£| — > 0, all systems 
are equivalent. This is to say 

Proposition 3.2. There exists positive constants S m ,(3 and C depending on 
the data a, b, c, d and v, such that for |£| < 5 m and for any t > 0, 

(3.2) || exp(-L4)|| < Cexp(-(3ft). 

Consequently, for any initial data Yq with supp(Yo) C [— S m ,d m ], 

E(t)<C t-^WYoWl,. 

Proof. By referring to the definitions of A, a and e, one sees that as 

iei-o, 

A 4£ 2 and tr(A) = va + e ~ {y + 1)£ 2 . 

Therefore, there exists 5 m > such that for £ in [— 5 m ,5 m ], A < and 

1 tr(A) 

- < < 2. 

2 - (u+i)e - 

()3.2j) then follows promptly from ()2.12j) . 

Now, for any initial data Yq with supp(lo) C [— S m , S m ], 

E(t) = [ \Y(t,^)\ 2 dC < C [ exp(-2/?i£ 2 K(sup(|^(0| 2 )) 

by using the change of variable r = \flj5i^. □ 

3.2. High frequency analysis. The complete dissipation and the partial 
dissipation cases have to be studied separately. In the latter case, we will 
give one example where the decay rate can be arbitrarily small. 

Introducing the number 



1, if r ^ 0, 
0, if r = 0, 



{r} 

for r£l. Then order(cr) = {a} + {c} - {6} - {d} and order(e) = 2 - 2{d}. 
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3.2.1. The complete dissipation case (v = 1). It is observed in the following 
that order(cr) dictates if the system is in the KdV-Burgers class or in the 
BBM-Burgers class. 

Proposition 3.3. Assume v = 1. For any 8 > 0, there exists (3 > ; such 
that if supp(Y ) C R\[-S, 8], 

(3.3) E{t) <exp(-2/?i) ll*o ||£| , 

for any t > 0. In addition, 

• if order(o~) < 0, the system is in the BBM-Burgers class. Namely, 
there exist positive constants 8m, (3 and C, such that for |£| > 8m 
and t > 0, 

||exp(-L4)|| < Ce -/W ; 

• if order(a) > 1, the system is in the KdV-Burgers class. Namely, 
there exist positive constants 8m, (3 and C, such that for |£| > 5m 
and t > 0, 

||exp(-L4)|| < Ce-^ H . 
Proof. ^From (|2.3j) . one finds that for £ almost everywhere, 

lj t \Y(t,0\ 2 + a(£)m,t)\ 2 + s(0\w(t,0\ 2 = 0. 

This gives directly, by setting (3 = mm{a(8),s(8)} which is positive, that 
1)3.3(1 is valid. Furthermore, || exp(— tA)\\ < Ce~^ 1 for |£| > 8. To figure 
out if the system is in the BBM-Burgers or in the KdV-Burgers class, we 
separate the cases as follows. 

• Assume first order(cx) > 1. Then either d = or b = 0. Without 
loss of generality, let us assume d = 0. 

— If A = (a — e) 2 — 4£ 2 <t 2 > for |£| large enough, then order(c) = 
1. In that case, there exist (3 > and 8m > 

det(A) _ a + a 2 a 

Al " W " JTT * m 

for |£| > 8 M - By using (l2~T31) 

|| exp(-L4)|| < C(l + t£ 2 )e- 2 ^ 2 < C7e"^ 2 

for |£| > 5m and the system is in the KdV-Burgers class. 

— If A < for |£| large enough, then there exists 5m > such 
that for |f| > 5 A /, < tr(A) < 2£ 2 , and (l2~T2|l implies that 
the system is in the KdV-Burgers class. 

• Assume now that order(cr) < 0. 

— If b 7^ and d ^ (weakly dispersive systems) then for |£| large 
enough, Re(Ai) < tr(A) ~ ^ + 4 as |£| — > oo. This shows that a 
damping like e~^ 2 is unlikely for high frequencies. Therefore 
the weakly dispersive systems are in the BBM-Burgers class. 
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if (6 ^ and d = 0) or (b = and d ^ 0). Without loss of 
generality, let us consider the case 6/0 and d = 0. Since 
A ~ £ 4 as |£| — ► oo, we have 

2det(A) 2(q + q 2 ) 

^(AT " -fTT ~ C " 

as |£| — > oo. This shows that a damping like e~^ 2 is again 
unlikely for high frequencies. Therefore the system is in the 
BBM-Burgers class. 

□ 



3.2.2. The partial dissipation case (y = 0). We first note that when a system 
satisfies (C2) hypothesis, a = and therefore Ai = at £ = a~2. But one 
can always chose 6m large enough so for |£| > <5m> o is positive, bounded 
from below and away from zero. Therefore the point where a vanishes will 
be considered in the next subsection. We now prove that in the partial 
dissipation case, the decay rate is related to order(cx) and the strength of 
the dissipation which is characterized by order (e). 

Proposition 3.4. With v = 0, 

• if order[a) > 2 — ^order(e), then the system is in the KdV-Burgers 
class; 

• if | order[a) \ < 2 — ^order(e), then the system is in the BBM-Burgers 
class. 

In above two cases, when Yq is supported in R\[— 6m, 6m], then for 
any t > 0, 

E(t)<Cexp(-2pt)\\Y \\ 2 Ll ; 

• if order(o~) < —2 + ^order(e), arbitrarily slow decay can occur. 

Proof. To begin, one observes that A = e 2 — 4£ 2 <r 2 . When A > 0, A! 
satisfies 

(3.4) 2Ai = £ I 1 - \ 1 - -^f- 

which is a direct consequence of (j2.8j) . 

• When order(e)=0, i.e d ^ (and order(a) < 1): 

— if order(cx) > —1, and if A = e 2 — 4£ 2 ct 2 > for |£| large enough, 
which is possible only for order(cr) = —1, we have 

det(A) 22 
Ai > -^y >Cia = 0(1) 

as |£| — > oo. Then by ()2.13j) we have 

\\e- tA \\ < C{l + t\i\a)e- tXl < Ce'^ 

for |£| large enough and the system is in the BBM-Burgers class. 
On the other hand, if A < for high frequencies, since tr(A) ~ 



12 



M. CHEN & O. GOUBET 



g as |£| — > oo, then (j2.12j) implies that the system is in the 
BBM-Burgers class; 

— if order(a) = -2, A ~ jp > as |£| -> oo and by (|3~1|) . 
Ai ~ C|£| -2 , therefore arbitrarily slow decay could occur. An 
example of such case will be given below. 

• When order(e) = 2 i.e d = (and order(a) > —I): A = £ 2 (£ 2 — 4cr 2 ) 
has a limit Ao in [— oo, +oo] when |£| approaches +oo. 

— If Ao is in (— oo,0], then since tr(A) = £ 2 , (|2,12|) implies that 
the system is in the KdV-Burgers class. This occurs when 
order(cr) = 2 and may occur when order(cr) = 1. 

— If Ao is in (0, +oo], then since 

2det(A) ^ ^ det(^) 2 

(|2.13j) implies for any £ 
(3.5) ||exp(-L4)|| < C[ 1 + |£|crmin (t, -U ) ) exp(-ta 2 ). 



A, . 

* If order (a) = 1, then (j3,5j) implies the system is in the 
KdV-Burgers class. 

* If order(<r) = 0, ^= = 0(1) as |£| — > oo, the system is in 
the BBM-Burgers class. And similarly, 

* if order (cr) = — 1, any arbitrarily slow decay could occur. 

□ 

Example of slow decay: Consider the linearized BBM-BBM system 
with partial dissipation, 

Vt + Ux ~ brixxt = 0, 
iM Vx du X xt — 
which has order(e) = and order(cr) = —2. Since as |£| —* +oo, 

A~I>0, 1,1=2151^-1=^ 

\ V (MA))'V / tr(A) d{2 



Therefore 

l|e--||<Cexp(-p) 

which shows that any arbitrary slow decay could occur. 

3.3. Middle range frequency analysis. We first note from Lemma 12.71 
that to get the optimal decay estimate for the cases where det (^4) (and there- 
fore Ai) has a zero for |£| > 0, these cases need to be discussed separately. 
Therefore, we have the following two propositions. 
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Proposition 3.5. Assume that v = 1, or that v = and the dispersive 
coefficients a,b,c,d satisfy (CI). Then for any 5 m and 5m, < 5 m < 5m, 
there exists ft > suc/i i/iat /or |£| G [6~ m , 5m] for any t > 0, 

(3.6) ||exp(-L4)|| < Cexp(-/3i). 

Moreover for QTiy Yq with suppovt included in [^m; ^iw] U [ — $M ? — ^m]? 

£7(i) < Cexp(-2/3t)||yo|li|. 

Proof. Since tr(A) and det(A) cannot vanish for |£| G [<5 m ,5jw] under the 
assumptions, ()3.6j) is the direct consequence of ()2.12l) and (|2.13j) . In addition 

E(t) < sup \\e- tA \\ 2 \\Y \\ 2 L2 <Cexp(-2/3i)||F ||| 2 , 

£ ? x 

which completes the proof of the proposition. □ 

Remark 3.6. By noticing that \3. b]) can be replaced by 

||exp(-L4)|| < Cexp(-/3*i£ 2 ). 

with f3* = (3/ 5 M , the middle range frequency analysis and the high frequency 
analysis can be combined to simplify certain calculations for these systems 
regardless if they are in BBM-Burgers class or KdV-Burgers class. 

Proposition 3.7. Assume that v = and the dispersive coefficients satisfy 
(C2). Then for any < 5 m < 5m with r = a~2 G [5 m ,5M], there exists 
(3 > and C > such that for any |£| G [6 m , 5m] and for any t > 0, 

(3.7) || exp(-L4)|| < Cexp{-/?t(|£| - r) 2 )}. 

Moreover for any Yq with support included in [5 m , 5m] U [— 5m , — 5 m ] and for 
any t > 0, 

E(t) < Ct-^WYoWli. 

Remark 3.8. Proposition \'J. 1\ shows that even when the dichotomy is not 
valid, the energy could decay as 0(t -1 / 4 ) when t goes to +oo. 



Proof. Since det(^4) vanishes at r = (y/a) 1 and, when |£| — > r, A ~ 
5^3 > 0, Ai ~ (3det{A) ~ pa 2 ~ | - r) 2 . Therefore, from (I7T31 . 

|| exp(— tA)\\ < C(l + min(i, l)o~)exp (-/?cr 2 t) 

as |£| in the neighborhood of r. Using the fact that for t > 0, min(t, 1) < 
so there exists C > such that 

|| exp(-L4)|| < C(l + \/to-)exp(-/5o- 2 t) < Cexp(-^o- 2 t), 

we obtain the estimate (|3.7j) for |£| close to r. For other |£| in [£ to ,<Sm]j the 
same argument in the proof of Proposition 13.51 and Remark 13 .61 applies . For 
the decay rate of E(t), same argument as in the proof of Propositions 13.21 
can be used. In fact, |£| — r plays the same role as |£| in that case. □ 
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3.4. Decay for linear systems. Since linear system (J1.3J) defines a semi- 
group e~ tA for t > 0, that is contracting on L 2 x L 2 in the variable (rj,w), 
the initial value problem is therefore well-posed and the L 2 norm decays. 

Combining the low, middle and high frequency analysis, the decay rate 
for the linear system ()1.3|) can be stated as 

Theorem 3.9. For systems hl.'J\) with the dispersive constants a, b, c, d sat- 
isfy the constraints (CO)-(Cl) or (C0)-(C2), assuming either {v = 1} or, 
{u = and order(a) > — 2 + ^order(s)} , then for any (t)o,Huq) = (r]o,wo) G 

D L 2 (M)) 2 where Huq = ( (i_"|2)(x+b| a ) )^o> there exists a constant 
C , such that for any t > 

\\( V ,Hu)\\ L 2 = \\(r!,w)\\ Ll <CrV\ 

Remark 3.10. This is equivalent to say , with respect to physical variables 
(rj,u), that for any (r] ,u ) £ [L 2 X n L l x ) x (H x n Wi 1 ' 1 ), 

\\{r},u)\\ LlxH u<Cr l / A 

for any t > 0, where h = order(H) = {a} + {d} — {c} — {&}. 

Proof. Combining the low, middle and high frequency analysis, we have 

E(t) = E low {t) + E middle (t) + E high {t) < C(r]o,u )(t~^ + e~ 2f3t ) 

for t > where C(rjo,uo) is a function of the dispersive coefficients and the 
norms of t/q and Uq. □ 
We complete this section by the following 

Corollary 3.11. For any dissipation, the classical Boussinesq system, the 
Bona-Smith system, the coupled KdV-BBM (b = c = 0) system, the BBM- 
KdV systems (a = d = 0) and the weakly dispersive systems (b > and 
d > ) with a < or c < belong to the BBM-Burgers class. 

Corollary 3.12. With complete dissipation, the KdV-KdV system (b = 
d = 0, a = c > 0) belongs to the KdV-Burgers class; the weakly dispersive 
systems (b > and d > 0) belong to the BBM-Burgers class. 

Remark 3.13. When the consideration is restricted to the linear systems, 
a result similar to Theorem W.fA (substituting a for e in the statement) holds 
when one replaces iyr\ xx ,u xx ) by ij] xx ,vu xx ) in the right-hand side of 
since r\ and u play the same role. 



4. Nonlinear Theory 

For convenience, we will only consider in this section (i) the complete 
dissipation and (ii) the partial dissipation with a,b,c,d satisfy the (CI) 
assumption. The partial dissipation with a, b, c, d satisfy (C2) will be con- 
sidered elsewhere. 
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4.1. A general result. Consider an evolution equation that reads 

(4.1) v t + Lv + d x (F(v)) = 0, v{t = 0) = v 

where L is a linear unbounded operator with symbol A and F is a nonlinear 
quadratic operator. 

Assuming that L generates a semi-group S(t) on l? x that satisfies the 
dichotomy assumption (|3.1j) . namely there exist (5 > and 5 > (5 can be 
+oo) such that for any t > and 

• for |e| < 5, \\S(t)\\ <Cexp(-/?t£ 2 ), 
•for |£|>5, ||S(t)|| <Cexp(-/ft). 

In addition, 

(4.3) sup^l^^n ) < dKIUinXg 

t>o 

Assuming also the nonlinear term F(v) satisfies 



/ \ 1/2 

(4.4) sup \F\+ ( / |£| 2 |F| 2 d£ I < G||« 

iei<<5 V-'i^ 5 / 



2 1 Z7l 2 I ^ /"Ml„. II 2 

L2 



for any t > 0, where -F = J-(F(v)). 

Let us recall that a mild solution to <|4.1jl is a solution to the integral 
equation 

(4.5) v(t) = S(t)v - [ S(t- s)d x F{v{s))ds. 

Jo 

Under the above assumptions, we may construct a solution to (|4.5|) by per- 
forming a fixed point argument (see ^Uj, [Tl], JJ) on the space 



E= lu: suplt 1 / 4 !!^)!!^} <oo), 



which is a Banach space of functions that are continuous in time with value 
in L 2 that are 0{t~ 1 ^) when t goes to +oo. If C is small enough, a fixed 
point argument to the Duhamel's form of the equation in the ball in E 
centered at origin would provide the solution. 

Theorem 4.1. For system j4-l\ ) with assumptions \4-ty - $T~3[ )- \4-4\ )> there 
exists a numerical constant C such that for any mild solution to \4-l\l start- 
ing from vq with 

then 

(4.6) \\v{t)\\ L 2 < 0{t- lfi ) as t^oo. 
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Proof. To begin, we first control the low frequency part of the nonlinear 
term. Let 

rt 

N:=T[ I S{t- s)d x F(v{s))ds 



Using the first inequality in (|4.'2[) in combination with (|4.4|) . one obtains 



\ 1/2 ( r 

/ \N\ 2 dA <C I [ \\e~^ A \\ 2 e\F\ 2 d( 
J\£\<6 J Jo [J\^\<S 



1/2 



ds 



(4.7) 



< C 



< c 



re 



1/2 



sup \F\ ds 

vl«l<* / 



t \\v{s)f 



% ds. 



We now control the high frequency part of the nonlinear term, using the 
second inequality in (|4.2I) in combination with (|4.4|) 

\ V2 „, / „ \ 1/2 



(4.8) 



\N\ 2 d£ 



< C 



\t\>6 



-0(t-s) 



f\F\ 2 d£ 



ds 



\t\>8 



< C / e-W- s )\\v(s)\\l 2 ds. 



Introducing the norm 
(4.9) 



M(t)= sup (s l ^\\v(s)\\ L 2), 

s£[0,t] 



and if v solves then due to (|4~7|) -- (j4~%|) . 

* 1/4 |K*)IU? <t 1/4 ||S(t>o||L ? +t 1/4 ||^)llL? 



<t 1/4 ||S(t)i;o|U2 +CM(t) 5 



+ _^exp(-/3(t- S )) 



s l/2( t _ g )3/4 s l/ 

By applying the change of variable s = rt in the integration, one finds 

t l/4 t l/4 



+ 



■ exp(-/?(t - s)) 



ds 



a l/2( t _ a )3/4 s l/2 

<C+ / — exp(-/3t(l -r))dr < C. 

Therefore, using the property (|4.3|) . the positive, nondecreasing function 
M(t) satisfies M(0) = and for any t > 0, 

(4.10) C M(t) 2 - M{t) + C > 0, 
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where Co is a positive constant. Choosing C such that Cqx 2 — x + C = 
has two real roots < n < ri, namely choosing C < then 1)4. lOj) holds 
only if M{t) is trapped in the interval [0, r{\. Therefore when 

n ,i ^ C 1 



Ci 4C C! 

M(i) is bounded and is valid. □ 

4.2. Applications to weakly dispersive systems with complete dis- 
sipation or partial dissipation. Since b > and d > 0, order(cr) < 
and the corresponding linearized system is in the BBM-Burgers class. From 
Proposition 13.31 and Proposition 13.41 this corresponds to consider complete 
dissipation, or partial dissipation together with order (<r) > — 1. 

Theorem 4.2. Consider a weakly dispersive two-way wave model (b > 
and d > 0) with either the complete dissipation or the partial dissipation 
together with a < or c < 0. Then, for small initial data, 

• if H is of order 0, 

(4-11) h{t)\\l % + \Ht)\\ll < Oir 1 / 2 ); 

• if H is of order 1, 

(4.i2) N(<)ll^ + K*)llfli<o(*" l/a ); 

• if H is of order —1, 
\m t i " 1 ' ' ! i l - 



(4.13) l|r?(t)|lL + ll^)HL<0(t- 1/2 ); 



as t — > oo. 

Proof. Note that the theorem is proved after (|4,4j) is validated and we 
will do that by discussing the cases according to the order of H. 

• If H is of order or 1. Introducing the change of variable 

v = Hu) = JF -I (r/, w), 

the full nonlinear system 

. , „, Vt + u x + au xxx - bn xxt + (rju) x = vr\ xx , 

(4.14) 

u t + rjx + c^a; - du xxt + uu x = u xx , 
transforms to 

v t + Lv = -d x F(v), 
where L has symbol A = ( . V ,° L N> ^ s S ri ( a - , i)? (J j an( ^ ^ reads 



F(v) 



(i - bdiy^H-^ 



•iv 



^Hil-dd^-^H^w) 2 



To check (|4.4|) . it is natural to separate the estimate into two parts. 
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— Low frequency < 5) estimate: Since H' 1 , which has order 
or —1, is bounded on L 2 , straightforward computations lead 
to 

WijfPhvW^ + iKH^yw^ < c([|/?[li s + 

<C(\\v\\h x + \\w\\ 2 L2 J. 

Because (1 — bd 2 ) -1 and H(l — dd 2 )~ l are bounded operators, 

sup \F\ < C\\vf L 2. 
\t\<8 

— High frequency > 5) estimate: Since d x (l — bd 2 )~ l is a 
smoothing operator (it is of order —1) and H -1 is a bounded 
operator on L 2 ., 

kcm^wh^wW^ < c(U\ 2 Ll + \\wWl,), 

where Lemma 2.2(ii) in pj, 

Wf9\\H--<C\\f\\ L 2\\g\\ L 2 

is used. Now, consider the term d x (H(l — dd%)~ ){(H w) 2 ) in 
F. If H is of order 0, it can be bounded in the same way. If H 
is of order 1, then 

\\d x H(l - ddlr^H-'wfW^ < C\\(H- l wf\\ Ll 

< C\\H- l wf Hl < C\\w\\ 2 Ll , 

where Lemma 2.2(iv) in [S] 

\\fg\\L0<C\\f\\ m \\g\\ m 

is used. 

Combining the lower frequency and higher frequency analysis, one 
sees IJ4.4JI is valid and Theorem 14. II yields the desired result. 
• If H is of order —1, introducing the change of variable 

(4.15) v = J 7 ~ 1 {(H- 1 ii,u)). 

and setting r = H~ l r], the full nonlinear system (|4.14j) reads 

vt + Lv = -d x F(v) 
where L has symbol A = ( . i^ct isgn(u!i)£,a\ 

The proof is then very similar to the previous case and therefore 
omitted. 

□ 
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Corollary 4.3. For the following two special cases, we have 

• solutions to Bona-Smith system (a = 0, 6 > 0, c < and d > 0) with 
complete or partial dissipation satisfy M.L'Jp ; 

• solutions to BBM-BBM system with complete dissipation (a = c = 
0,1/ = 1) satisfy UU\) . 

4.3. Application to KdV-Burgers systems with complete dissipa- 
tion. Using Proposition 13.31 this implies that order(<r) > 1, and then that 
6 = and/or d = 0. 

First case: Consider the case where 6 = d = 0. The analysis in jS] implies 
that a = c = 1/6, so the system satisfies (C2) assumptions. Since the di- 
chotomy assumption Q4.2j) was proved in Section |3] and the linearized system 
is in the KdV-Burgers class, Theorem 14.11 applies when (|4.4[) with 5 = +oo 
is verified. 

Let us observe that H = 1 and that the full nonlinear system reads 

vt + Lv = -d x F(v) 

where v = (v,u), L has symbol A = ( . f ... i s g n -(ui)£,°~ \ an( ^ 

m = (| 

Since 

\\F(v)\\Ll<C(\\v\\ 2 Ll + \\u\\ 2 Li ), 
(|4.4f) with 5 = oo is a direct consequence. 

Second case: Consider the case where b > and d = 0. Due to (CO), 
c > 0. Since order(cr) > 1, the system must have a = c > 0. Therefore 
order {H) = — 1. Introducing the change of variable (|4.15f) the system reads 
as (j4.3|) with the following nonlinearity 

, . /H- 1 (l-bd 2 xy 1 (uHT) 

F {V) = 

V 2 

(|4.4|) with 5 = oo is valid and the proof is straightforward and then omitted. 

Third case: Consider the case where 6 = and d > 0. Since order (c) = 1, 
then a and c can not vanish and order(i?) = 1. In this case, with the change 
of variable ()4.2|) . the system reads as (|4.3|) with F being 



iHil-ddZy^H^w) 2 
(|4.4(1 with 5 = oo is again valid and the proof is straightforward and then 
omitted. 

Therefore we can state 



Theorem 4.4. Consider a KdV-Burgers system with complete dissipation. 
Then for small initial data, 



20 



M. CHEN & O. GOUBET 



• if H is of order 0, 

(4.16) n#)iiii + imi)iii|<or 1/2 ); 

• if H is of order 1, 

(4.17) umk + Ht)\\ 2 Hi < o(t-^ 2 y, 

• if H is of order —1, 

(4.18) Ht)\\ 2 m x + \\u(t)\\ll<0{t~ l l 2 ); 
as t — > oo. 

Remark 4.5. This is not surprising for KdV-KdV system (b = d = 0) 
since if we introduce the new variables rj = w 1 + w 2 and u = w 1 — w 2 , 
then \4-14\ ) reads as a system of two linear KdV-Burgers systems (weakly) 
coupled through nonlinear terms. See Section 2.3 in jS]. 

4.4. Other cases. In some other cases, the method presented here does not 
work straightforwardly. As pointed out in the Introduction, other methods 
exist which might enable the analysis to go further. These methods could 
also be helpful to extend our local results to global ones. 



5. The L£°-decay rate 

First, we observe that for the cases of weakly dispersive wave equations 
and KdV-KdV system, the nonlinear terms satisfy 

/ \ 1/2 

(5.1) sup|£F|+(/ \^\F\ 2 dA <C\\v\\ l1 \\v x \\lI- 

We now estimate the decay rate of d x v{t) in L 2 when v solves ()4.5(l . To 
begin, we differentiate (|4.5|) with respect to x and treat the nonlinear term 
of the resulting equation with a procedure similar to the one in the proof of 
Theorem 14.11 but using ()5.1|) instead of ()4.4|) . 

We first note that for low frequencies, 

V/iei<<5 / Jo (i-s) J / 4 

and for high frequencies 

( f e\N\ 2 dt) <C [\-W- s ^v\\ L2 Jv x \\ Ll ds. 

\J\t\>6 J Jo 



We now consider the linear part. For vq in H 1 (R) n L 2 (R) , the linear term 
can be estimated by splitting the region of integration into low frequencies 
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and high frequencies. Using (|4.2j) . we have 
\\d x S(t)v \\ 2 L 2 < C 
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ee-^)\\v \\h+e- 



<r n- II 2 _i_ r'c- 2 / 3 *!!,, II 2 



XW01II2 



Therefore, the linear part behaves like 0(t 3//4 ) as t — > 00. Since [|v ||^2 is 
0(i -1 / 4 ) as t — > 00, we have 

ll^WIki <cWr 3 / 4 + c snp (s l ^\\v(s)\\ L2 ) 

se[o,t] v 7 

t e -/3(t-s) 







Simple calculations show 

'■ l ds 



L 



+ 



S l/4(t- s )3/4 y Q s l/4 



ds < C. 



(5.2) 



Therefore 

||^(t)|| L 2 <C(v )t~ 3 ^ + C sup ||^( S )|| Li sup (^^(^ii^) 

se[o,t] se[o,t] 

<C(v )r 3 ^ + C 2 M(t) sup ||9^(s)|| L 2. 

ss[0,t] 

Since M(t), that is defined in (|4.9|) . is bounded by the first root n of 

C x 2 -x + C = 0, 

and n ~ C as C -» (since ^(1 - (1 - 4C C) 1/2 ) ~ 2C as C -» 0). 
Therefore, when C small enough, M(t) < 2C. Hence by choosing C such 
that 



2C 2 C < ~ 



(5.3) 
l|5~2J) leads to 

(5.4) \\d x v(t)\\L* <2C(v )t- 3 / 4 

and we obtain the following theorem. 

Theorem 5.1. For system \4-l\) with assumptions i4-%!\ )- n~3[ )- \4-4\ )> assume 
Vq is in H 1 (M) nL 1 (M) and ||fo||L 1 nL 2 ^s small enough. Then 

NU« < o(^ 1/2 ) 

as t — > 00. 

Proof. Using (j5,4[) and (j4.6j) together with 

I/2 II^Wiixv 



II / II ll 1 / 2 !!^ Ill/ 2 

v\\L2> < IMI/2 \\o x v\ f 



yields the desired result. 



□ 
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6. Numerical Result 

Numerical simulations are performed on several systems and results on 
BBM-BBM and Bona-Smith systems with complete or partial dissipations 
are reported here. The results show not only that the theoretical results on 
the decay rates are sharp, but also the constants involved are reasonably 
sized. 

In these numerical computations, the initial data are taken to be 

T] = sech 2 (^(x - xo)), 
u = Vo- ?7o/ 4 > 

where xq is in the spatial domain [0, L], where L is taken to be large enough 
so the solution near the boundary is smaller than the machine roundoff error 
during the whole computation. The spectral method is used on the spatial 
domain [0, L\ and the leap-frog algorithm is used on the time advancing. 
The decay rate r and the constant C in 

\\v || ~ Ct~ T , as t — > oo 

is calculated by first computing 

lo Sj7 

t-n— 1 

The computation is stopped when r{t n ) is approaching to a constant and 
the value r is obtained by averaging the last 5 data. The constant C is then 
computed by averaging the last five |M|(i n )i^. 

In the computations reported below, L = 320, dx = 0.1 and dt = 0.05, 
where dx and dt are the meshsize in space and time respectively. 

BBM-BBM system (a = c = 0,b = d= 1/6) with complete dissipa- 
tion. It is shown in Theorem 14.21 and 15. II that for small data, 



< Ci^ 1/4 and |H| L oc < C 2 t- 1/2 . 



The numerical computation is performed for time interval [0,50], and the 
result shows 

|M| L2 ~ 1.4232t- a2470 and |M| L °c ~ 1.4989^ a4963 . 

Therefore, it is clear that the theoretical result is sharp and the constants 
involved are not large. Moreover, it seems that the small data requirement 
might be removed if, for example, other methods were employed. 

Bona-Smith system (a = 0, b = — c = d = 1/3) with complete and 
partial dissipation. This case is again covered by Theorem 14.21 and 15.11 

By direct computation, we obtain for complete dissipation, 

\\v\\ L2 ~ 1.4015t- a2477 , \\v\\ L oo ~ 1.4466t- a4998 , 
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and for partial dissipation 

1Mb ~ 0.6676t- a2519 , || V || £ =o ~ 0.6595t- a5105 . 
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